Energy scales and the non-Fermi liquid behavior in YbRh2Si2 
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Multiple energy scales are detected in measurements of the thermodynamic and transport prop- 
erties in heavy fermion metals. We demonstrate that the experimental data on the energy scales 
can be well described by the scaling behavior of the effective mass at the fermion condensation 
quantum phase transition, and show that the dependence of the effective mass on temperature and 
applied magnetic fields gives rise to the non-Fermi liquid behavior. Our analysis is placed in the 
context of recent salient experimental results. Our calculations of the non-Fermi liquid behavior, 
of the scales and thermodynamic and transport properties are in good agreement with the heat ca- 
pacity, magnetization, longitudinal magnetoresistance and magnetic entropy obtained in remarkable 
measurements on the heavy fermion metal YbRh2Si2. 

PACS numbers: 71.27.-l-a, 71.10.Hf, 73.43.Qt 



An explanation of the rich and striking behavior of 
strongly correlated electron ensemble in heavy fermion 
(HF) metals in the vicinity of a quantum phase transi- 
tion is, as years before, among the main problems of the 
condensed matter physics. It is common wisdom that 
low-temperature and quantum fluctuations at quantum 
phase transitions form the specific heat, magnetization, 
magnetoresistance etc. which are drastically different 
from those of ordinary metals [H-Q • Conventional argu- 
ments that quasiparticles in strongly correlated Fermi liq- 
uids "get heavy and die" at the QCP commonly employ 
the well-known formula basing on assumptions that the 
z-factor (the quasiparticle weight in the single-particle 
state) vanishes at the points of second-order phase tran- 
sitions 0. However, it has been shown this scenario is 
problematic Q. On the other hand, facts collected on 
HF metals demonstrate that the effective mass strongly 
depends on temperature T, doping (or the number den- 
sity) X and applied magnetic fields B, while the effective 
mass M* itself can reach very high values or even di- 
verge, see e.g. d, Such a behavior is so unusual that 
the traditional Landau quasiparticles paradigm does not 
apply to it. The paradigm says that elementary excita- 
tions determine the physics at low temperatures. These 
behave as Fermi quasiparticles and have a certain effec- 
tive mass M* which is independent of T, x, and B and 
is a parameter of the theory Q. 

A concept of fermion condensation quantum phase 
transition (FCQPT) preserving quasiparticles and inti- 
mately related to the unlimited growth of A/*, had been 
suggested [gt-fllj. Studies show that it is capable to de- 



liver an adequate theoretical explanation of vast majority 
of experimental results in different HF metals [l^, [3 • 
contrast to the Landau paradigm based on the assump- 
tion that M* is a constant, in FCQPT approach M* 
strongly depends on T, x, B etc. Therefore, in accord 
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with numerous experimental facts the extended quasi- 
particles paradigm is to be introduced. The main point 
here is that the well-defined quasiparticles determine as 
before the thermodynamic and transport properties of 
strongly correlated Fermi-systems, M* becomes a func- 
tion of T, X, B, while the dependence of the effective 
mass on T, x. B gives rise to the non-Fermi liquid (NFL) 
behavior ^U^. 

In this letter, we analyze the NFL behavior of strongly 
correlated Fermi systems and show that this is gener- 
ated by the dependence of the effective mass on tem- 
perature, number density and magnetic field at FCQPT. 
We demonstrate that the NFL behavior observed in the 
transport and thermodynamic properties of HF metals 
can be described in terms of the scaling behavior of 
the normalized effective mass. This allows us to con- 
struct the scaled thermodynamic and transport proper- 
ties extracted from experimental facts in wide range of 
the variation of scaled variable. We show that "pecu- 
liar points" of the normalized effective mass give rise to 
the energy scales observed in the thermodynamic and 
transport properties of HF metals. Our calculations of 
the thermodynamic and transport properties are in good 
agreement with the heat capacity, magnetization, lon- 
gitudinal magnetoresistance and magnetic entropy ob- 
tained in remarkable measurements on the heavy fermion 
metal YbRh2Si2 For YbRh2Si2 the constructed 

thermodynamic and transport functions extracted from 
experimental facts show the scaling over three decades in 
the variable. 

To avoid difficulties associated with the anisotropy 
generated by the crystal lattice of solids, we study the 
universal behavior of heavy-fermion metals using the 
model of the homogeneous heavy-electron (fermion) liq- 
uid [ii,|2i!,|2|. 

We start with visualizing the main properties of FC- 
QPT. To this end, consider the density functional theory 
for superconductors (SCDFT) 23]. SCDFT states that 
at fixed temperature T the thermodynamic potential $ is 
a universal functional of the number density n(r) and the 
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anomalous density (or the order parameter) K(r,ri) and 
provides a variational principle to determine the densi- 
ties [i^] ■ At the superconducting transition temperature 
Tc a superconducting state undergoes the second order 
phase transition. Our goal now is to construct a quantum 
phase transition which evolves from the superconducting 
one. In that case, the superconducting state takes place 
at T = while at finite temperatures there is a normal 
state. This means that in this state the anomalous den- 
sity is finite while the superconducting gap vanishes. For 
the sake of simplicity, we consider a homogeneous Fermi 
(electron) system. Then, the thermodynamic potential 
reduces to the ground state energy E which turns out 
to be a functional of the occupation number n{p) since 
K = y/ n{l — n) [TgI [23l - [25j . Upon minimizing E with 
respect to n(jp), we obtain 
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where /i is the chemical potential. It is seen from Eq. ([T]) 
that instead of the Fermi step, we have < n{p) < 1 in 
certain range of momenta pi < p < Pf with k is finite in 
this range. Thus, the step-like Fermi filling inevitably un- 
dergoes restructuring and formes the fermion condensate 
(FC) as soon as Eq. ([T]) possesses not-trivial solutions at 
some point x = Xc when Pi = Pf = Pf [E, 12, 13]. Here 
Pi? is the Fermi momentum and x = p^/Sn'^. 

At any small but finite temperature the anomalous 
density k (or the order parameter) decays and this state 
undergoes the first order phase transition and converts 
into a normal state characterized by the thermodynamic 
potential $o- At T 0, the entropy S = -d^o/dT of 
the normal state is given by the well-known relation Q 

5o = -2 1 [n(p) ln(n(p)) + {1 - n(p) ln(l - n{p))]^, 

which follows from combinatorial reasoning. Since the 
entropy of the superconducting ground state is zero, it 
follows from Eq. ^ that the entropy is discontinuous at 
the phase transition point, with its discontinuity AS" — 
Sq. The latent heat q of transition from the asymmetrical 
to the symmetrical phase is q — TcSq — since Tc = 
0. Because of the stability condition at the point of the 
first order phase transition, we have $o[n-(p)] = $[k(p)]. 
Obviously the condition is satisfied since q = 0. 

At T = 0, a quantum phase transition is driven by a 
nonthermal control parameter, e.g. the number density 
X. To clarify the role of x, consider the effective mass 
M* which is related to the bare electron mass m by the 
well-known Landau equation which is valid when M* 
strongly depends on _B, T or a; [l6j 
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Here we omit the spin indices for simplicity, n(p, T) is 
quasiparticle occupation number, and F is the Landau 



amplitude. At T = 0, Eq. ([3]) reads [H, [27 
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Here A^o is the density of states of free electron gas and 
F^{x) is the p-wave component of Landau interaction 
amplitude F. When at some critical point x = Xc, F^{x) 
achieves certain threshold value, the denominator in Eq. 
(HI) tends to zero so that the effective mass diverges at 
T = [ii-ii]. It follows from Eq. gj that beyond the 
quantum critical point (QCP) Xc^ the effective mass be- 
comes negative. To avoid unstable and physically mean- 
ingless state with a negative effective mass, the system 
must undergo a quantum phase transition at QCP x = Xc 
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FIG. 1: Schematic phase diagram of the system driven to 
the FC state. The number density x is taken as the control 
parameter and depicted as x/xc- The quantum critical point 
(QCP), x/xc = 1, of FCQPT is shown by the arrow. At 
x/xc < 1 and sufficiently low temperatures, the system is in 
the Landau Fermi liquid (LFL) state as shown by the shadow 
area. At T = and beyond QCP, x/xc > 1, the system is 
at the quantum critical line depicted by the dash line and 
shown by the vertical arrow. The critical line is characterized 
by the FC state with finite superconducting order parameter 
K. At Tc — 0, the order parameter k is destroyed, the system 
undergoes the first order phase transition and exhibits the 
NFL behavior at T > 0. 

Schematic phase diagram of the system which is driven 
to FC by variation of x is reported in Fig. [1] Upon ap- 
proaching the critical density Xc the system remains in 
LFL region at sufficiently low temperatures [l^, [l^l : that 
is shown by the shadow area. At QCP Xc shown by the 
arrow in Fig. [l] the system demonstrates the NFL be- 
havior down to the lowest temperatures. Beyond QCP 
at finite temperatures the behavior is remaining the NFL 
one and is determined by the temperature-independent 
entropy So [Si- In that case at T — 0, the system is ap- 
proaching a quantum critical line (shown by the vertical 



arrow and the dashed hne in Fig. [T]) rather than a quan- 
tum critical point. Upon reaching the quantum critical 
line from the above at T ^> the system undergoes the 
first order quantum phase transition, which is FCQPT 
taking place at Tc = 0. 

At T > the NFL state above the critical line, see 
Fig. [U is strongly degenerated, therefore it is captured 
by the other states such as superconducting (for exarnple, 
by the superconducting state in CeCoIns l2ll [22l . [23 |) 
or by AF state (e.g. AF one in YbRh9Si9 [l6j|l lifting 
the degeneracy. The application of magnetic field B > 
Bco restores the LFL behavior, where Bco is a critical 
magnetic field, such that at i? > Bco the system is driven 
towards its Landau Fermi liquid (LFL) regime [2^ . In 
some cases, for example in HF metal CeRu2Si2, -B^n = 0, 
see e.g. while in YbRhaSia, BcO - 0.06 T (So]. In 

our simple model Bco is taken as a parameter. 

Schematic phase diagram of the HF metal YbRh2Si2 
is reported in Fig. [2] Magnetic field B is taken as the 
control parameter. The FC state and the region lying 
at x/xc > 1, see Fig. (TJ can be captured by the super- 
conducting, ferromagnetic, antiferromagnetic (AF) etc. 
states lifting the degeneracy Since we consider 

the HF metal YbRh2Si2 the AF state takes place 
as shown in Fig. [5J As seen from Fig. [51 at elevated 
temperatures and fixed magnetic field the NFL regime 
occurs, while rising B again drives the system from NFL 
region to LFL one. Below we consider the transition re- 
gion when at rising B the system moves from NFL regime 
to LFL one along the dash-dot horizontal arrow, and at 
elevated T it moves from LFL regime to NFL one along 
the solid vertical arrow. 

To explore a scaling behavior of M*, we write the 
quasiparticle distribution function as ni(p) = n{p,T) — 
n(p), with n(p) is the step function, and Eq. ([3|) then 
becomes 
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At QCP the effective mass M* diverges and Ec^. ([5]) 
becomes homogeneous determining M* as a function of 
temperature 



M*(T) oc 



(6) 



while the system exhibits the NFL behavior [1^ |l5| . If 
the system is located before QCP, M* is finite, at low 
temperatures the system demonstrates the LFL behav- 
ior that is M*(T) ~ M*+aiT'^, with oi is a constant, see 
the inset to Fig. [2] Obviously, the LFL behavior takes 
place when the second term on the right hand side of Eq. 
([5]) is small in comparison with the first one. Then, at ris- 
ing temperatures the system enters the transition regime: 
M* grows, reaching its maximum M^j a,t T — Tm, with 
subsequent diminishing. Near temperatures T > Tm the 
last "traces" of LFL regime disappear, the second term 
starts to dominate, and again Eq. ^ becomes homoge- 
neous, and the NFL behavior restores, manifesting itself 
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FIG. 2: Schematic phase diagram of the HF metal YbRh2Si2. 
Bco is magnetic field at which the effective mass diverges. 
AF denotes antiferromagnetic (AF) state. At B < B^o the 
system is in AF state. The vertical arrow shows the transition 
from the LFL regime to the NFL one at fixed B along T 
with M* depending on T. The dash-dot horizontal arrow 
illustrates the system moving from NFL regime to LFL one 
along B at fixed T. The inset shows a schematic plot of 
the scaling behavior of the normalized effective mass versus 
the normalized temperature. Transition regime, where 
reaches its maximum value M^j at T = Tm, is shown by the 
hatched area both in the main panel and in the inset. The 
arrows mark the position of inflection point in and the 
transition region. 



in decreasing M* as T ^Z'^. When the system is near 
QCP, it turns out that the solution of Eq. @ M*{T) can 
be well approximated by a simple universal interpolating 
function [l^; d; 113 ■ The interpolation occurs between 
the LFL (M* ~ M* + a^T^) and NFL (M* oc T'^/S) 
regimes thus describing the above crossover [l^, [iBl • In- 
troducing the dimensionless variable y = Tn = T /Tm, 
we obtain the desired expression 
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Here = M* /Ml^ is the normalized effective mass, 
Co = (1 -I- C2)/(l -I- ci), ci and C2 are fitting parameters, 
parameterizing the Landau amplitude. 

The inset to Fig. [2] demonstrates the scaling behavior 
of the normalized effective mass M'^ = M* /M^j ver- 
sus normalized temperature T^r = T /Tm, where M^,j is 
the maximum value that M* reaches at T = Tm- At 
T < Tm the LFL regime takes place. At T > Tm 
the r^2/3 j-ggijjig takes place. This is marked as NFL 
one since the effective mass depends strongly on tem- 
perature. The temperature region T ~ Tm signifies the 
transition between the LFL regime with almost constant 
effective mass and NFL behavior, given by T-2/3 de- 
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pendence. Thus temperatures T ^ Tm can be regarded 
as the transition region between LFL and NFL regimes. 
The transition temperatures are not reahy a phase tran- 
sition. These necessarily are broad, very much depending 
on the criteria for determination of the point of such a 
transition, as it is seen from the inset to Fig. [2] As 
usually, the transition temperature is extracted from the 
temperature dependence of charge transport, for exam- 
ple, from the resistivity p{T) — + AT^ with po is the 
residual resistivity and A is the LFL coefficient. The 
crossover takes place at temperatures where the resis- 
tance starts to deviate from the LFL behavior. Ob- 
viously, the measure of the deviation from the LFL 
behavior cannot be defined unambiguously. Therefore, 
different measures produce different results. 

It is possible to transport Eg. Q to the case of the 
application of magnetic fields [H, EE 113 ■ The applica- 
tion of magnetic field restores the LFL behavior so that 
depends on B as 

Mil « - B,o)-^/\ (8) 

while 

TmxMB-B,o), (9) 

where /is is the Bohr magneton [l^ [l^, [l^. Employ- 
ing Eqs. dH) and dU to calculate Mlj and Tm, we 
conclude that Eq. ^ is valid to describe the normal- 
ized effective mass in external fixed magnetic fields with 
y = T/{B - Bco)- On the other hand, Eq. © is valid 
when the applied magnetic field becomes a variable, while 
temperature is fixed T = Tf. In that case, as seen from 
Eqs. ([6]), (O and®, it is convenient to rewrite both the 
variable as y = (_B — Bco)/Tf, and Eq. © as 

Pb{Bm ~ B,o) xTf. (10) 

It follows from Eq. (O that in contrast to the Landau 
paradigm of quasiparticles the effective mass strongly de- 
pends on T and B. As we will see it is this dependence 
that forms the NFL behavior. It follows also from Eq. ([7]) 
that a scaling behavior of M* near QCP is determined 
by the absence of appropriate external physical scales to 
measure the effective mass and temperature. At fixed 
magnetic fields, the characteristic scales of temperature 
and of the function M* (T, B) are defined by both Tm and 
respectively. At fixed temperatures, the character- 
istic scales are {Bm — Bcq) and M^.j. It follows from Eqs. 
(|H]) and (0) that at fixed magnetic fields, Tm — > 0, and 
M^^ — > CO, and the width of the transition region shrinks 
to zero as i? — !> Bcq when these are measured in the exter- 
nal scales. In the same way, it follows from Eqs. ([6]) and 
([TU| that at fixed temperatures, {Bm — Bco) ~^ 0, and 
M^j — ^ oo, and the width of the transition region shrinks 
to zero as ^0. Thus, the application of the external 
scales obscure the scaling behavior of the effective mass 
and thermodynamic and transport properties. 

A few remarks are in order here. As we shall see, 
magnetic field dependencies of the effective mass or of 




FIG. 3: The normalized entropy Sn{B / Binf) versus y = 
B/Binj and the normalized entropy Sn{T /Tinf) versus y = 
T /Tinf calculated at fixed temperature and magnetic field, 
correspondingly, are represented by the solid lines and shown 
by the arrows. The inflection point is depicted by the dash- 
dot arrow. 



other observable like the longitudinal magnetoresistance 
do not have "peculiar points" like maximum. The nor- 
malization are to be performed in the other points like 
the inflection point at T = T^f (or at B = Binf) shown 
in the inset to Fig. [Hby the arrow. Such a normalization 
is possible since it is established on the internal scales, 
Tinf oc Tm oc (_B - Sco). 

In what follows, we compute the effective mass and 
employ Eq. ([7]) for estimations of considered values. To 
compute the effective mass M*{T,B), we solve Eq. ([5]) 
with special form of Landau interaction amplitude, see 
Refs. [H, llEl for details. Choice of the amplitude is 
dictated by the fact that the system has to be at QCP, 
which means that first two p-derivatives of the single- 
particle spectrum £(p) should equal zero. Since first 
derivative is proportional to the reciprocal quasiparti- 
cle effective mass 1/M*, its zero just signifies QCP of 
FCQPT. Zeros of two subsequent derivatives mean that 
the spectrum e(p) has an inflection point at pp so that 
the lowest term of its Taylor expansion is proportional 
to (p - pf)^ [ill- After solution of Eq. (O, the ob- 
tained spectrum had been used to calculate the entropy 
S{B,T), which, in turn, had been recalculated to the 
effective mass AI* (T, B) by virtue of well-known LFL re- 
lation M*{T,B) = S{T,B)/T. Our calculations of the 
normalized entropy as a function of the normalized mag- 
netic fleld B / Binf = y and as a function of the normal- 
ized temperature y = T/Tinf are reported in Fig. [3) Here 
Tinf and Binf are the corresponding inflection points in 
function S. We normalize the entropy by its value at 
the inflection point SN{y) = S{y)/S{l). As seen form 
Fig. [21 our calculations corroborate the scaling behavior 
of the normalized entropy, that is the curves at different 
temperatures and magnetic fields merge into single one 
in terms of the variable y. The inflection point Tinf in 
S{T) makes M*{T, B) have its maximum as a function of 
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T, while M*{T, B) versus B has no maximum. We note 
that our calculations of the entropy confirm the validity 
of Eq. ([7|) and the scaling behavior of the normalized 
effective mass. 

Exciting measurements of C/T oc AI* on samples of 
the new generation of YbRh2Si2 in different magnetic 
fields i? up to 1.5 T [l^ allow us to identify the scaling 
behavior of the effective mass M* and observe the dif- 
ferent regimes of M* behavior such as the LFL regime, 
transition region from LFL to NFL regimes, and the NFL 
regime itself. A maximum structure in C'/T oc M^^ at 
Tm appears under the application of magnetic field B 
and Tm shifts to higher T as i? is increased. The value 
of C/T = 7o is saturated towards lower temperatures 
decreasing at elevated magnetic field, where 70 is the 
Sommerfeld coefficient [l8| . 

The transition region corresponds to the temperatures 
where the vertical arrow in the main panel of Fig. [2] 
crosses the hatched area. The width of the region, being 
proportional to Tm cc {B — Bco) shrinks, Tm moves to 
zero temperature and 70 ex. M* increases as B ^ Bco- 
These observations are in accord with the facts llSll. 
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by the hatched area in the inset to Fig. [2l are clearly seen 
in Fig. m illuminating good agreement between the the- 
ory and facts. It is worthy of note that the normalization 
procedure allows us to construct the scaled function C /T 
extracted from the facts in wide range variation of the 
normalized temperature. Indeed, it integrates measure- 
ments of C /T taken at the application of different mag- 
netic fields into unique function which demonstrates the 
scaling behavior over three decades in normalized tem- 
perature as seen from Fig. 2] 
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FIG. 5: The field dependencies of the normalized magneti- 
zation M collected at different temperatures shown at right 
bottom corner are extracted from measurements collected on 
YbRu2Si2 19]. The kink (shown by the arrow) is clearly seen 
at the normalized field Bn ~ B/Bk — 1. The soUd curve 
represents our calculations. 

Consider now the magnetization M as a function of 
magnetic field B at fixed temperature T — Tf 



M{B,T)= / x{b,T)db, 



(11) 



FIG. 4: The normalized effective mass M'^ extracted from 
the measurements of the specific heat C/T on YbRh2Si2 in 
magnetic fields B 18] hsted in the legend. Our calculations 
are depicted by the solid curve tracing the scaling behavior of 



To obtain the normalized effective mass M'^ , the max- 
imum structure in C/T was used to normalize C/T, and 
T was normalized by Tm- In Fig. |4]M^ as a function of 
normalized temperature T/v is shown by geometrical fig- 
ures, our calculations are shown by the solid line. Figure 
U] reveals the scaling behavior of the normalized exper- 
imental curves - the scaled curves at different magnetic 
fields B merge into a single one in terms of the normalized 
variable y = T/Tm- As seen, the normalized mass 
extracted from the measurements is not a constant, as 
would be for LFL. The two regimes (the LFL regime and 
NFL one) separated by the transition region, as depicted 



where the magnetic susceptibility x is given by 

I3M*{B,T) 



X{B,T) = 



1 



(12) 



Here, /3 is a constant and Fq is the Landau amplitude 
related to the exchange interaction Q. In the case of 
strongly correlated systems > -0.9 (Hd^j. There- 
fore, as seen from Eq. ([T^ . due to the normalization the 
coefficients f3 and (1 + Tq ) drops out from the result, and 
X oc M*. 

One might suppose that Fq can strongly depend on B. 
This is not the case, since the Kadowaki- Woods ratio is 
conserved 30 32], A{B)/-f^{B) oc A{B)/x'^iB) oc const, 
we have 70 oc M* oc x- Here A is the coefficient in the 
T^ dependence of resistivity p. 

Our calculations show that the magnetization exhibits 
a kink at some magnetic field B = Bk- The experimental 
magnetization demonstrates the same behavior [loj . We 
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use Bk and M{Bk) to normalize B and M respectively. 
The normalized magnetization M{B)/M{Bk) extracted 
from facts [l^ depicted by the geometrical figures and 
calculated magnetization shown by the solid line are re- 
ported in Fig. [5j As seen, the scaled data at different 
Tf merge into a single one in terms of the normalized 
variable y = B/T^. It is also seen, that these exhibit 
energy scales separated by kink at the normalized mag- 
netic field Bpf = B/Bk — 1. The kink is a crossover 
point from the fast to slow growth of M at rising mag- 
netic field. It is seen from Fig. [SJ that our calculations 
are in good agreement with the facts, and all the data 
exhibit the kink (shown by the arrow) at Bat ~ I taking 
place as soon as the system enters the transition region 
corresponding to the magnetic fields where the horizon- 
tal dash-dot arrow in the main panel of Fig. [5] crosses 
the hatched area. Indeed, as seen from Fig. [5l at lower 
magnetic fields M is a linear function of B since M* is 
approximately independent of B. Then, it follows from 
Eqs. d?]) and ^ that at elevated magnetic fields M* be- 
comes a diminishing function of B and generates the kink 
in M(B) separating the energy scales discovered in Rcfs. 
[13, [3| • Then, as seen from Eq. ([TU)) the magnetic field 
Bk at which the kink appears, Bk — Bm '^Tf, shifts to 
lower B as TV is decreased. This observation is in accord 
with facts [lElll. 
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FIG. 6: Magnetic field dependence of the normalized mag- 
netoresistance pM versus normalized magnetic field, pjv was 
extracted from LMR of YbRh2Si2 at different temperatures 
p^l. [l9| listed in the legend. The infiection point is shown by 
the arrow, and the soHd line represents our calculations. 

Consider a longitudinal magnetoresistance (LMR) 
p[B,T) = po + AT"^ as a function of B at fixed Tf. 
In that case, the classical contribution to LMR due to 
orbital motion of carriers induced by the Lorentz force 
is small, while the Kadowaki- Woods relation [30l - l3^ , 
K = A/7g oc A/y^ = const, allows us to employ M* 
to construct the coefficient A [sl], since 70 oc x M* . 
As a resuh, p{B,T) - po oc (M*)^. Fig. [6] reports the 
normalized magnetoresistance 



versus normalized magnetic field y = B/Binf at differ- 
ent temperatures, shown in the legend. Here pinf and 
Binf are LMR and magnetic field respectively taken at 
the inflection point marked by the arrow in Fig. jS] Both 
theoretical (shown by the solid line) and experimental 
(marked by the geometrical figures) curves have been nor- 
malized by their inflection points, which also reveals the 
scaling behavior - the scaled curves at different temper- 
atures merge into single one as a function of the variable 
y and show the scaling behavior over three decades in 
the normalized magnetic field. The transition region at 
which LMR starts to decrease is shown in the inset to 
Fig. |2]by the hatched area. Obviously, as seen from Eq. 
P0|) . the width of the transition region being propor- 
tional to Bm — Binf oc Tf decreases as the temperature 
Tf is lowered. In the same way, the infiection point of 
LMR, generated by the inflection point of AI* shown in 
the inset to Fig. [5] by the arrow, shifts to lower B as 
Tf is decreased. All these observations are in excellent 
agreement with the facts [O, • 

The evolution of the derivative of magnetic entropy 
dS{B,T)/dB as a function of magnetic field B at fixed 
temperature Tf is of great importance since it allows us 
to study the scaling behavior of the derivative of the ef- 
fective mass TdM*{B,T)/dB oc dS{B,T)/dB. While 
the scaling properties of the effective mass M* {B, T) can 
be analyzed via LMR, see Fig. jB] 




PNiy) 



pjy) - Po 

Pinf 



(13) 



123456789 
NORMALIZED MAGNETIC FIELD 

FIG. 7: Normalized magnetization difference divided by tem- 
perature increment {AM/AT)n versus normalized magnetic 
field at fixed temperatures listed in the legend is extracted 
from the facts collected on YbRh2Si2 [131 • Our calculation 
of the normalized derivative [dS/dB)N — {AM / AT)n versus 
normalized magnetic field is shown by the solid line. 

As seen from from Eqs. Q and (|T0| . at y < 1 the 
derivative —dMN{y)/dy oc y with y = (B — Bco) / {Bmf — 
Bco) oc {B — Bco)/Tf. We note that the effective mass as 
a function of B does not have the maximum. At elevated 
y the derivative —dMjs[{y) / dy possesses a maximum at 
the inflection point and then becomes a diminishing func- 
tion of y. Upon using the variable y = {B — BcQ)/Tf, 
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we conclude that at decreasing temperatures, the lead- 
ing edge of the function —dS/dB oc —TdM*/dB becomes 
steeper and its maximum at {Binf — Bco) oc Tf is higher. 
These observations are in quantitative agreement with 
striking measurements of the magnetization difference di- 
vided by temperature increment, — AAZ/AT, as a func- 
tion of magnetic field at fixed temperatures Tf collected 
on YbRh2Si2 t20|]. We note that according to the well- 
know thermodynamic equality dM/dT = dS/dB, and 
AM /AT ~ dS/dB. To carry out a quantitative analy- 
sis of the scaling behavior of —dM*{B,T)/dB, we cal- 
culate the normalized entropy S shown in Fig. |3] as a 
function of B/Binf at fixed temperature Tf. Fig. [Tire- 
ports the normalized (dS/ dB) jv as a function of the nor- 
malized magnetic field. The scaled function {dS/dB)^ 
is obtained by normalizing {—dS/dB) by its maximum 
taking place at Bm, and the field B is scaled by Bm- 
The measurements of —AM /AT are normalized in the 
same way and depicted in Fig. [7] as {AM/AT)^ versus 
normalized field. It is seen from Fig. [7] that our calcula- 
tions shown by the solid line are in good agreement with 
the facts and the scaled functions {AM/ AT)m extracted 
from the facts show the scaling behavior in wide range 
variation of the normalized magnetic field B/Bm- 




B(T) 

FIG. 8: Temperature versus magnetic field T — B phase di- 
agram for YbRh2Si2. Solid circles represent the boundary 
between AF and NFL states. The solid squares denote the 
boundary of the NFL and LFL regime [l7l.Il9l. [sH shown by 
the dotted line which is approximated by \/B — Bco Di- 
amonds mark the maximums Tm of C/T [3l shown in Fig. 
m The dash-dot line is approximated by Tm oc a{B — Bco), a 
is a fitting parameter, see Eq. Q. Triangles along the solid 
line denote Ti„f in LMR \VA [l9l | sown in Fig. [T] the solid 
line represents the function Tinf oc b{B — Bco), & is a fitting 
parameter, see Eq. pup . 



Fig. [8] reports T^f and Tm versus B depicted by the 
solid and dash-dotted lines, respectively. The bound- 
ary between the NFL and LFL regimes is shown by the 
dashed line, and AF marks the antiferromagnetic state. 
The corresponding data are taken from Ref. (l7l - [T9l . [3ll | . 
It is seen that our calculations are in good agreement 
with the facts. In Fig. [51 the solid and dash-dotted lines 
corresponding to the functions T^f and Tm, respectively, 
represent the positions of the kinks separa ting the energy 
scales in C and M reported in Ref. [13, Il9l | It is seen 
that our calculations are in accord with facts, and we 
conclude that the energy scales are reproduced by Eqs. 
dH) and ([T0| and related to the peculiar points T^f and 
T]\i of the normalized effective mass M^ which are shown 
by the arrows in the inset to Fig. [^j 

At B -> BcQ both Tinf ^ and Tm ^ 0, thus the 
LFL and the transition regimes of both C/T and M 
as well as these of LMR and the magnetic entropy are 
shifted to very low temperatures. Therefore due to ex- 
perimental difficulties these regimes cannot be often ob- 
served in experiments on HF metals. As it is seen from 
Figs. [31[5l[ni[7|and[Sl the normalization allows us to con- 
struct the unique scaled thermodynamic and transport 
functions extracted from the experimental facts in wide 
range of the variation of the scaled variable y. As seen 
from the mentioned Figures, the constructed normalized 
thermodynamic and transport functions show the scaling 
behavior over three decades in the normalized variable. 

In summary, we have analyzed the non-Fermi liquid 
behavior of the heavy fermion metals, and showed that 
extended quasiparticles paradigm is strongly valid, while 
the dependence of the effective mass on temperature, 
number density and applied magnetic fields gives rise to 
the NFL behavior. We have demonstrated that our the- 
oretical study of the heat capacity, magnetization, longi- 
tudinal magnetoresistance and magnetic entropy are in 
good agreement with the outstanding recent facts col- 
lected on the HF metal YbRh2Si2. Our normalization 
procedure has allowed us to construct the scaled thermo- 
dynamic and transport properties in wide range of the 
variation of the scaled variable. For YbRh2Si2 the con- 
structed thermodynamic and transport functions show 
the scaling behavior over three decades in the normal- 
ized variable. The energy scales in these functions are 
also explained. 
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